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Abstract 

We analyse the effects of large-scale inhomogeneities upon the observables of a grav- 
itational lens system, focusing on the issue of whether large-scale structure imperils the 
program to determine the Hubble parameter through measurements of time delays be- 
tween multiple images in lens systems. We find that the lens equation in a spatially 
flat Robertson- Walker cosmology with scalar metric fluctuations is equivalent to that 
for the same lensing system in the absence of fluctuations, but with a different angular 
position of the source relative to the lens axis. Since the absolute position of the source 
is not observable, gravitational lens measurements cannot directly reveal the presence 
of large scale structure. Large-scale perturbations do not modify the functional rela- 
tionship between observable lens parameters and the Hubble parameter, and therefore 
do not seriously affect the determination of Hq from lens time delays. 
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1 Introduction 

Gravitational lenses are proving to be both intrinsically fascinating systems and valuable 
astrophysical tools that can help determine the fundamental cosmological parameters (for a 
review see Schneider et al. 1992). In principle, lens observations can provide estimates of the 
Hubble parameter H Q and the density parameter Q Q (Refsdal 1964). The consistency between 
gravitational lens estimates and more local determinations of these parameters would also 
constitute additional evidence in favor of the standard Friedmann-Robertson- Walker (FRW) 
cosmological models. 

Measurements of time delays between multiple images in a gravitational lens system 
were first shown by Refsdal (1964) to provide a potentially direct determination of H Q , 
independent of the standard cosmic distance ladder. If the structure of the deflector in a 
gravitational lens system is sufficiently well understood and its redshift known, then the 
relationship between the redshifts of the images and the deflector, and the images' relative 
angular positions, relative magnifications, and time delays, provide a determination of the 
deflector parameters, such as its mass, as well as of the Hubble parameter. To date, the 
double quasar 0957+561, the first lens system discovered, is the only one for which a time 
delay has been reliably measured (Vanderriest et al. 1989, Lehar et al. 1992, Press et al. 
1992), and a value for H derived (Falco et al. 1991, Rhee 1991, Roberts et al. 1991). For this 
system, the estimated errors in H Q are around 10% to 30%, due to uncertainties in the density 
parameter fl Q , in the relative angular positions of the images and the deflecting galaxy, and 
especially in the lens model parameters. It is reasonable to expect that measurements of 
time delays in other gravitational lens systems that have simpler structure than 0957+561 
will improve this estimate in the near future, and provide a valuable tool for cosmology. 

The aim of this article is to discuss whether the existence of large-scale structure in the 
Universe compromises the program to use time delay measurements to determine the Hubble 
parameter. In principle, inhomogeneities in the mass distribution along the photon paths 
affect the observables of a gravitational lens system, and if not properly taken into account 
as part of the model for the lens system, they could systematically bias the determination of 
H Q from lens time delays. Among others, Alcock & Anderson 1985, 1986, Watanabe, Sasaki, 
& Tomita 1992, Sasaki 1993, and Seljak 1994 have estimated the effect of departures from a 
FRW cosmology on the time delay between multiple images in a lens system. (Other recent 
studies of light propagation in perturbed FRW cosmologies include Durrer 1994 and Pyne 
and Birkinshaw 1994.) 

There are two effects of large-scale structure on the relation between lens time delays 
and H Q . First, long-range fluctuations in the gravitational potential near the line of sight 
to a lensed QSO can affect the distance measure, and thus the proportionality between the 
measured lens-induced time delay and H Q . The time delay is proportional to D^Dg/D^ ~ 
H' 1 , where D d and D s are the distances from the observer to the lens deflector and the 
source respectively, and D ds is the distance between deflector and source. Here, the central 
issue of debate is what distance D is appropriate: the FRW angular-diameter distance, the 
Dyer-Roeder (Dyer and Roeder 1972) distance (for an empty or partially filled beam), or 
some other distance measure which also takes into account the effects of large-scale shear. 
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The second effect of large-scale perturbations is a direct contribution to the lens time delay, 
in addition to that arising from the lens itself. (This direct contribution may be present even 
if the lens-induced delay vanishes in the absence of perturbations; it is therefore conceptually 
distinct from the indirect distance- measure effect.) 

This paper is primarily concerned with the latter, direct effect of perturbations on lens 
time delays. Our main conclusion is that, while large-scale inhomogeneities do have an 
effect upon time delays and other observables in a gravitational lens system, they do not 
compromise the program to extract from them the value of H Q . More precisely, the lens 
equation in the presence of scalar metric fluctuations is the same as that describing an 
identical lens system in the absence of fluctuations but with an (unobservably) different 
absolute angular position of the source. It is thus observationally impossible to distinguish 
time delays induced by large-scale structure from intrinsic delays due to the lens itself. The 
important corollary is that, to leading order in the fluctuations, the relationships between 
the observables of the lens system, the Hubble parameter, and the lens model parameters in 
the presence of scalar metric fluctuations are the same as if the inhomogeneities were absent, 
modulo distance measure effects. We also show that for large-scale, small-amplitude density 
perturbations, the modification of the distance measure from the FRW angular-diameter 
distance is small. 

Taken together, these results imply that large-scale structure does not imperil the pro- 
gram to determine H Q from time delay measurements in gravitational lenses. On the other 
hand, it also implies that one cannot use time delay measurements to detect or constrain 
large-scale inhomogeneities in the Universe, once H Q is determined reliably by other means. 
It was suggested a few years ago (Allen 1989) that time delay measurements in gravitational 
lenses could serve as gravitational wave detectors. The same technique could in principle 
have been extended to probe large-scale inhomogeneities in the matter distribution (Frieman 
& Turner 1989, unpublished). We have argued, however, (Frieman, Harari, & Surpi 1994) 
that it is not observationally possible to distinguish the time delay induced by gravity waves 
(tensor metric fluctuations) from the intrinsic time delay originating in the lens geometry. 
Here we extend this result to scalar metric fluctuations of a FRW spacetime (arising from 
matter-density fluctuations), using the same technique based upon Fermat's principle in 
curved space-time. Our conclusions are similar to, but our methods differ from, those of 
Frieman, Kaiser, & Turner 1990 (unpublished). Our main result is that the lens equation 
in the presence of scalar metric fluctuations is effectively the same (to the relevant order of 
approximation) as in the absence of fluctuations, but with the lens system having a different 
alignment between observer, deflector, and source. 



2 The Lens Equation and Large-scale Structure 

Consider a thin, stationary gravitational lens, embedded in a spatially-flat FRW cosmology, 
with scalar metric perturbations representing large-scale matter-density fluctuations. We 
assume a weak gravitational field U(x,t) for the deflector in the lens system and small- 
amplitude metric inhomogeneities, and we work to first order in both the large-scale metric 
fluctuation amplitude h^ u and the deflector potential U (which implies first order in a, 
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the deflection angle imprinted on the light rays). We use the longitudinal or conformal- 
Newtonian gauge, and let Greek indices fi, v run from to 3 while Latin indices i,j run from 
1 to 3; we also set the speed of light c = 1. Defining the conformal time drj = dt/a(t), with 
t the proper time and a(t) the cosmic scale-factor, the metric is 

ds 2 = a 2 (r])[(l + 2{U + <f)))dr] 2 - (1 - 2{U + cj)))5 ij dx i dx j ] , (1) 

where <p(x, if) completely describes the scalar metric fluctuation for a matter-dominated 
universe, and is equal to a gauge-invariant variable; it satisfies the relativistic generalization 
of the Poisson equation. For Q — 1, the conformal time and scale-factor can be normalized so 
that t = 2n 3 /3H Q and a{rf) = 2r] 2 /H , where r] — 1 and subscript zero indicates the present 
epoch. For the growing mode of adiabatic scalar density fluctuations in a matter-dominated 
(0 = 1) universe, the metric perturbation amplitude is time-independent, (f>(r],x) = 4> a {x) 
(e.g., Mukhanov et al. 1992). The results below can easily be extended to the non-flat case 

i- 

Without loss of generality, we can place the observer at the origin of coordinates and the 
z- axis coincident with the lens axis (the line that joins observer and deflector); the deflector 
is located at D = (0, 0, r d ) and the source at S = (x s , 0, r s ), see Figure 1. To avoid confusion 
with redshift, we use r to denote the value of the coordinate z along the lens axis. We 
denote angular positions by two-component vectors, defined by the (x, y) components along 
the photon paths in the plane perpendicular to the lens axis. The indices a, b, running from 
1 to 2, denote these components. Thus we characterize the absolute angular position of the 
source relative to the lens axis by the vector f3 = (x s /r s , 0) (these are the components at the 
source position, z = r s ). 

From Eq. ([I]), if the spatial photon trajectories are known, one can evaluate the conformal 
time of travel by integration along r (the coordinate along the lens axis), 



rr s 

7] Pd I dr 
Jo 



(2) 



The first three terms in the integrand of Eq. (0) are the geometric contributions to the travel 
time, while the last contains the gravitational potential contributions from the lens deflector 
and large-scale metric perturbations. 

We now determine the photon paths by implementing Fermat's principle (Blandford & 
Narayan 1986, Kovner 1990, Nityananda and Samuel 1992). We first determine the time of 
travel along null trial paths of the metric ([!]). Each trial path is composed of two segments, 
one from the source to the deflector plane, and another from the deflector plane to the 
observer, as appropriate for a thin lens. Each segment is a solution of the geodesic equations 
for the metric (|l|) neglecting the gravitational potential U of the deflector (which is taken 
into account through the bending at the deflector), but including the effect of the large-scale 
metric fluctuation 0. We further require that the paths built in this way are null paths of the 
full metric ([[]) (i.e., they satisfy ds 2 = and Eq. (|2|)) but not necessarily geodesies. Along 
them, we evaluate the conformal travel time from source to observer. Fermat's principle 
states that null geodesies are those null paths for which the arrival time is an extremum. 
Making the conformal travel time an extremum leads to the lens equation, from which the 
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relationships among the apparent angular positions of the multiple images, their time delays, 
relative magnifications, and the lens parameters can be read off. 

The affine parameter for the photon geodesies of the metric (jl|), neglecting the deflector 
gravitational potential U, may be written as A = r + Q{4>). To leading order in <ft, the 
corresponding geodesic equations may then be integrated to obtain 

£ = « + A(r) , (3) 

(recall that x(r) denotes the vector whose (x,y) components give the photon trajectory 
parametrized in terms of the distance along the lens axis) where e denotes an arbitrary 
integration constant. For the scalar growing mode, A is given by 

A» = -2 f dr<f>, a , (4) 

where a comma denotes an ordinary derivative in the FRW metric. Here, we have assumed 
that these trajectories form small angles with the lens axis, and we work to first order in 
and in the trajectory angle from the lens axis. That is, we can expand the function in eqn.(|4]) 
around the lens axis, A a (x b ,r) = A a (r)+A a b x b (r) + ..., where A a (r) = A a (0,r) is the function 
evaluated along the lens axis. Since A a (r) oc a , the approximation in eqns.^) and (|), i.e., 
dropping terms of order A a b , corresponds to keeping only first derivatives (gradients) of the 
potential, neglecting terms of order d 2 <p/ dx a dx h . This means that we are not including 
the relative focusing due to large-scale fluctuations. For typical gravitational lenses, the 
maximum transverse separation between the image paths is £ ~ DO ~ 10 — 20 kpc, much 
smaller than the wavelengths A ^ 10 Mpc of the large-scale perturbations we are interested 
in. Thus, the second derivative terms we are neglecting are suppressed compared to the 
first derivative terms by a factor £/A ^ 10~ 3 . We shall use this approximation consistently 
throughout. As a result, the integrand in Eq. (f§) is to be evaluated along the lens axis, 
which one can think of as a fiducial or 'mean' photon path. 

We now use Eq.([|), with appropriate values for the integration constants e in each seg- 
ment < r < r d and < r < r s , to build a zig-zag trajectory that starts at S, is deflected 
at r = rd (the deflector plane), and arrives at the observer at the origin. There is a family 
of trajectories satisfying these focusing conditions, which we choose to parametrize in terms 
of the apparent angular position 6 of the source image relative to the deflector (since this 
is an observable quantity). To evaluate 6, we use again Eq. ([|) with appropriate boundary 
conditions to evaluate the photon trajectory x d (r) that arrives directly from the deflector 
to the observer simultaneously with the source image. The apparent angular position of the 
image relative to the deflector as seen by the observer is then given by 

dx d 



dr 



(5) 

r=0 



r=o dr 

The family of trajectories that meet the focusing conditions at the source and the observer, 
parametrized in terms of 6, satisfies 

dx n ri A(r) , 

— = 6 - I — — dr + A(r) if r < r d 

dr Jo rd 

^ = -I*-0 + l± P- P—dr + A(r) if r > r d (6) 

dr r ds r ds Jr d r ds 
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where ra s = r s — r<±. The deflection angle imprinted by the lens upon the trajectory at the 
deflector plane (r = ra), which we denote by a, is given by 



a 



dx 

dr 



dx 

dr 



(7) 



Using Eq. (0), we find 



a = ^(0-/3 cff ) , 

rds 



where the effective misalignment angle (3 eS has been defined as 

f3 eS = f3 + f3 LSS , 



with 



(3 



LSS 



1 /" r d 1 

— / A(r)dr 



r d Jo 



r s Jo 



A.(r)dr 



(8) 



(9) 



(10) 



Written this way, Eqn. (§) is identical in form to the equation relating the image angular 
position with the source misalignment angle (3 and the deflection a in the absence of 
density fluctuations, but now in terms of an effective misalignment (3 eG . The geometric 
meaning of f3 eS is apparent from Eqn. (^): it is the source-deflector misalignment angle that 
an observer would infer from lens observations, assuming a homogeneous FRW spacetime, 
i.e., with no knowledge of large-scale perturbations. Alternatively, it is the source angular 
position that the observer would measure in the perturbed FRW universe in the limit that 
the lens mass vanishes (a — > 0). 

Now we evaluate the conformal time of travel by integrating Eq. (fj) along the null 
trajectories of Eq. @. It is useful to distinguish two contributions to the travel time, the 
geometric, i] g , and the potential, r] p . The former is easily evaluated: 



dr 



2\dr 2\dr 



r s + ~ r -^0- (0-2(3) 

2 r ds 



:n) 



Here, we have discarded a ^-independent term, irrelevant to extremizing the travel time. 
Note that Eq. (|TTJ) has the same form as the geometric contribution to the time of travel 
in the absence of density fluctuations. The fluctuations are taken into account, however, 
through the fact that 6 is the apparent angular position of the source image relative to the 
deflector in the presence of the metric perturbation 0. This is a crucial step that allows us 
to establish the equivalence between a lens in the presence of density fluctuations and a lens 
in the absence of fluctuations but with a different source alignment. 

The evaluation of the gravitational potential contribution to the time of travel requires 
more care, so we divide it into two parts, the first due to the large-scale metric fluctuations, 
r/p SS , and the second due to the potential U of the deflector, rj^ L . For the large-scale structure 
contribution, to first order in 6 and 9 we have 



^ P LSS 



-2 / dr<f)(x a ,r) 
Jo 



(12) 
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and the integration is taken along the unperturbed path, i.e., the trajectories of Eq. (^]) with 
A a = 0. In the integrand, we can expand the metric amplitude as 

<f)(x a , r) = 0(0, r) + 0, a x a (r) + <j), a , b x a (r)x b (r) + ... , (13) 

where x a (r) is given by integrating @. In keeping with the approximation used in eqns.(|3] 
and we neglect terms beyond the first derivative on the RHS of eqn.flTB]). Substituting 
(|j~3|) into fll2"D , integrating by parts, and retaining only # a -dependent terms, we find 

lss = _ e .( f d dr A - P dr r -±A) = -™ • (3 LSS , (14) 
\Jo Jr d r s J r ds 

with f3 LSS as defined in Eq. (0). As above, the middle terms in Eq. (|1"4]) are to be integrated 
along the fiducial mean path. 

The last contribution to the time of travel, due to the gravitational potential U of the 
deflector, is a local effect in a thin lens. It results in a function where £ is the impact 

parameter. Integrating eqn.(|5]), one finds £ = x(r d ) = r d 0. Thus, rj can be expressed in 
terms of the lens mass density projected on the lens plane, £(£) (Schneider, Ehlers & Falco 
1992), 



-2 £ U dr = -AG J d 2 i' £(£') In 



-m ■ (is) 



Putting together the results from Eqs. flnp, ( JT1 ) and (JI5[), the total conformal travel 

LSS _|_ /viG- 
P ~ 'P 



time for a source image, rj = r/ g + ?7p SS + r]^ L , becomes (up to ^-independent terms) 



V = r s + (0 - 2(3 cS ) ■ 6 - ij(r d 0) , (16) 



with f3 eS as defined in Eq. (pup. The main conclusion of this article derives from equation 



(|i6|), which may be thought of as an "equivalence principle" for gravitational lenses. The 
travel time (|T6"D is equivalent to the time of travel in a lens system identical to the one 



considered from the outset, but with no metric perturbations and with a different source 
alignment given by (3 eS instead of (3. However, as noted above, f3 eS is precisely the mis- 
alignment angle the observer would infer from lens observations without taking into account 
the large-scale perturbations. Thus, the observer can self-consistently ignore the presence of 
metric perturbations from the outset and obtain the correct time delay. The second term 
in equation (|16"D, which partially originates in the metric fluctuation 0, masquerades as an 
intrinsic geometric effect of the lens. The lens equations, deduced from the condition that 
rj be an extremum under variations of 9 for fixed source position, are the same in the two 
cases, so an observer measuring redshifts, relative angular positions, relative magnifications, 
and time delays, is unable to tell the two situations apart-we will show this explicitly below. 

To apply the results above, we must translate the conformal time delay between two 
images at angular positions #i and 2 , 

Arj = v{ei ) - 77(02) = ^(fc - 2 ) ■ (6, + 6 2 - 2/3 cff ) - (v(6) - V>(6)) , (17) 

into the proper time delay At measured by the observer. We treat the first and second terms 
on the r.h.s. of Eq. (|17D separately. The second term, due to the deflector gravitational 
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potential, is a purely local effect, which occurs when the photon paths are close to the 
deflector. Thus, the passage from conformal time-delay to proper time-delay at the observer's 
position is simply given by the factor (1 + z d ), with z d the deflector redshift, which accounts 
for the time dilation caused by the expansion of the Universe as the photons travel from the 
deflector to the observer. The first term in Eq. ([H]), on the other hand, is translated to the 
observer's proper time by multiplying it by the present scale factor a Q (a good approximation 
since At << H~ l ). In order to make contact with observations, it is best to express the 
result in terms of the source and deflector redshifts z s , z d . For the leading order contribution 
to the time delay, we evaluate the relationship between comoving coordinates and redshifts 
in the unperturbed FRW background (see discussion at the end of this section). For a photon 
emitted at coordinate distance r e at time t = t e and observed at r = at time t = t D in a 
matter-dominated spatially-flat FRW cosmology (e.g., Schneider et al. 1992), 

2 vm-i d. 



H a y/\ + z a e 



where the last term defines the angular-diameter distance D e to the event with comoving 
coordinate r e . In terms of angular-diameter distances, the observer's proper time delay 
between the two images becomes 



At = (1 + z d ) 



D>lDs '0i - e 2 ) ■ (^±^ - (3 cS ) - - ^(D d o 2 )) 



D 



ds 



(19) 



The appearance of the angular diameter distances D d , -Ds and D ds in this expression does 
not imply the need to independently determine the distance scale to the source and deflector: 
the combination (1 + z d )D d D B /D dB should simply be taken as shorthand for 

n ^ ,D d D s 2 (l-v/rT^Xl-y/TT^) 

( + ZA) D ds ~ H D VT+lT d - VT+Y S ■ (2U) 

The lens equation, obtained from the requirement that the time of travel for each image be 
an extremum with respect to variations in 9 (Fermat's principle) is now 

^ DlD '(0- M (21) 



do a 



ds 



This agrees with Eq. (|8|), since the deflection angle imprinted by the lens is given as a function 
of the lens mass distribution by a = dip/ 8$,. 

In the limit that the relative change induced by large-scale structure upon the angular 
separation between multiple images is small (|/3 LSS | << \0\ — 0% |), use of the lens equation 
(pl|) allows the time delay to be written in a simpler form, 

At = At intrinsic + At LSS , (22) 

where At mtrmsic is the time delay due to the lens geometry evaluated in the absence of metric 
perturbations, and At LSS is the lowest-order effect of the metric perturbations, 

At LSS = -(l + ^ d )%^/W(0i-02) • (23) 



Here it suffices to evaluate the angular separation between the images in the absence of 
metric perturbations. 



In the general case, we can use the lens equation (21) to rewrite the proper time delay 



for the image observed at angular position Oi in terms of lens observables as 



(1 + z d )D d 2D ds ddi 

where x(Q) = 4>{6)/Dd is the dimensionless lens potential. The fact that this expression does 
not contain (3lss, and is therefore independent of 0, embodies our main conclusion: while it 
is clear, for instance from Eqs. (|22|) and (|23|) , that large-scale inhomogeneities do influence 
the time delay between multiple images in a gravitational lens, at the same time they affect 
the other lens observables in such a way that they leave no observable tracks of their presence. 
That is, in Eq. ([24]), lens observations in principle provide the required distance factors and 
the angle 0i, while modelling of the lens (based on an inferred velocity dispersion and the 
observed surface brightness distribution) provides a parametrized determination of \i with 
no reference to large-scale structure. Thus, (|23|) is just the time delay that an observer with 
no knowledge of large-scale structure would use to constrain his or her lens model. 

The final thread in the argument that large-scale fluctuations do not perturb the relation 
between H Q and measured lens time delays involves the distance measure. In eq. (|T8|) and 
following, we used the angular-diameter distance for the FRW background. Since the matter 
in the universe along the line of sight to a lensed system is clumped, it is not clear that 
this is the appropriate distance measure to use, and there is a large literature that treats 
this thorny issue. Here, we note that by focusing on the effects of small-amplitude, large- 
scale perturbations, we may essentially sidestep this debate. For linear perturbations, one 
can expand the distance measure around the FRW background (e.g., Sasaki 1993). The 
distance measure perturbation is then proportional to the density perturbation averaged 
over the beam. It is clear that this average is generally less than of order the typical 
perturbation amplitude. Thus, for linear density perturbations, 5p/p -C 1, the distance 
measure perturbation is small, 5D/D ^ Sp/p <C 1. 



3 Conclusion 

The measurement of time delays between multiple images in a gravitational lens system can 
yield an estimate of the Hubble constant, as well as the lens mass or other lens parameters. 
Indeed, once a model for the deflecting object is assumed, knowledge of the redshifts of the 
images and the deflector, of their apparent relative angular positions, relative magnifications, 
and time delay allow a determination of the Hubble constant, since At cx DdD s /Dd s oc H~ l . 

For the lens equation in the presence of large-scale matter-density fluctuations, the per- 
turbations appear only through the quantity /3lss °f Eq. fllPf ), which is simply added to the 
intrinsic lens misalignment (3, to give the effective misalignment angle f3 cS = f3 + f3 LSS . Thus, 
matter-density fluctuations have exactly the same observational consequences as a change 



9 



in the (unobservable) misalignment angle between the source and the lens axis. The stan- 
dard method to determine the Hubble constant H Q proceeds from the lens observables in the 
usual way (Schneider et al. 1992), as if the large-scale inhomogeneities were absent. Indeed, 
the unobservable effective misalignment /3 eff cancels out from the expressions that relate the 
angular separation between the images, 6\ — 6 2 , and the time delay At (Cf. eqn. (plD ). 

Thus, large-scale matter-density fluctuations along the line of sight do not compromise the 
program to determine the value of the Hubble constant through time-delay measurements 
between multiple images in gravitational lenses. This also implies that lens time delay 
measurements are not likely to provide information about large-scale inhomogeneities in the 
matter distribution. 

It is important to emphasize that this conclusion is limited to the effects of small- 
amplitude (i.e., linear), large-scale density fluctuations. The assumption of linearity enters in 
several places. First, it implied that the scalar potential is time- independent in a spatially 
flat FRW background. While this simplifies the calculation, it is not an essential ingredient. 
More important, it is only for small-amplitude density perturbations that one can justify 
using the angular-diameter distance rather than a perturbed distance measure. This does 
not preclude the possibility that non-linear structure on small scales (e.g. galaxies) could 
significantly perturb the distance measure. The assumption of 'large scale' entered when we 
neglected second derivatives of the potential: this restricts the treatment to perturbation 
wavelengths A much larger than the maximum transverse path separation, £ ~ 10 — 20 kpc. 
Observations of galaxy clustering indicate that fluctuations in the galaxy density are non- 
linear on scales less than several Mpc, so the restriction to linearity already implies that the 
assumption A ^> £ is satisfied. Finally, we note that our conclusions are valid not only for 
scalar metric perturbations, but also for tensor metric fluctuations (gravitational waves), as 
shown in Frieman, Harari & Surpi 1994. 
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Figure 1: Gravitational lens geometry. The observer is at the origin of coordinates. The 
z-axis coincides with the lens axis, and the ^-coordinate is labeled by r to avoid confusion 

— * — * — * 

with redshift. The source forms an angle (3 with respect to the lens axis. Q\ and 9 2 are 
the apparent angular positions of the images and a is the deflection angle. 



